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A unified model is constructed to study the recently observed DNA entropic elasticity, cooperative 
extensibility, and supercoiling property. With the introduction of a new structural parameter (the 
folding angle <p), bending deformations of sugar-phosphate backbones, steric effects of nucleotide 
basepairs, and short-range basestacking interactions are considered. The comprehensive agreement 
of theoretical results with experimental observations on both torsionally relaxed and negatively su- 
percoiled DNAs strongly indicates that, basestacking interactions, although short-ranged in nature, 
dominate the elasticity of DNA and hence are of vital biological significance. 
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Recent in vitro experiments done on single double- 
stranded DNA (dsDNA) molecules or DNA-protein com- 
plexes reveal that DNA double-helix has nontrivial elas- 
tical properties |l|-|7j. At low external forces, it can be 
viewed as a simple wormlike chain Q ; at moderate forces, 
it becomes a rod with a large stretch modulus. But if 
pulled with a large force of about 70 pN the molecule 
as a whole can suddenly be driven to an almost fully 
stretched state with a contour length 1.7 times its native 
value HPJ . More strikingly, if a dsDNA at the same time 
has a slight deficit in linking number, i.e., negatively su- 
percoiled, a pulling force as small as 0.3 pN can distort 
the native structure of DNA considerablely ■ 

On the theoretical side, to understand these novel 
properties of DNA is of current interest and many mod- 
els are suggested and some valuable insights attained 
|p|,pij|-|20| . For example, to interpret the extensibility of 
DNA, some authors suggest one-dimensional two-state 
models with ||,[l2 13 1{| or without nearest-neighbor 
interactions; and to explain the supercoiling property of 
DNA, wormlike rod chain models, with (lj|l(| or with- 
out M,0,M bend-twist and/or stretch-twist coupling, 
are investigated. Nevertheless, a unified description still 
seems to be lacking and the underlying mechanism which 
should account for DNA cooperative extensibility and 
novel supercoiling properties still elusive. Here, we show 
that it is possible to understand all these experimental 
observations from a unified viewpoint. 

A simple elastical model is proposed by taking into ac- 
count the structural properties of realistic dsDNA. Bend- 
ing energy of the sugar-phosphate backbones, basestack- 
ing interaction between adjacent nucleotide basepairs as 
well as their steric effects on DNA axial bending rigidity 
are considered. We have introduced a new structural pa- 
rameter, the folding angle p. Model calculations on the 



elastical properties of torsionally relaxed and negatively 
supercoiled DNAs are in quantitative agreement with all 
the known experimental observations [||^| ■ The model 
indicates that basestacking interaction is the main fac- 
tor determining the highly extensibility and unwinding 
instability of DNA. We suggest that the present model, 
after some revisions, will also be able to account for the 
elasticity of positively supercoiled DNAs . 

In the model, the two inextensible backbones of DNA 
pl[ are characterized by the same bending rigidity k = 
ksT£ p , where £ p ~ 1.5 nm is their bending persis- 
tence length (BPL). Their position vectors are r, = 
J ti(s')ds', where t; (i — 1,2) is the unit tangential 
vector of the ith backbone, and s its arc length. The nu- 
cleotide basepairs between the backbones pl| are viewed 
as rigid planar structures with finite area and volume. 
Firstly, we consider bending energy of the backbones 
alone and each basepair connecting the two backbones 
is replaced for the moment by a thin rigid rod of length 
2R, with a unit vector b pointing along it from ri to 
r 2 , i.e., r 2 (s) — ri(s) = 2Rb(s). Relative sliding of 
the backbones is prohibited, the basepair planes are as- 
sumed to lie perpendicular to DNA central axis and 
b ti = b t 2 = |p2f. The central axis of dsDNA can be 
defined as r(s) = ri(s) + i?b(s), and its tangential vector 
is denoted by t, with t • b = 0. Since both ti and t 2 
lie on the same plane perpendicular to b, we obtain that 
ti = t cos p + n sin p and t 2 = tcost^ — rising, where 
n = b x t and p> is half the rotational angle from t 2 to 
ti (b being the rotational axis). We call <p the folding 
angle, it is in the range between — ir/2 and +7r/2 (p > 
for right-handed rotations and < for left-handed ones). 
It is not difficult to verify that 

db/ds= (t 2 -t x )/2R = -nsinp/R, (1) 

(here and after ds always denotes arc length element of 
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the backbones). With Eq. (|l|) and the definition of r we 
know that 



dv/ds = (ti + t 2 )/2 = tcosc/?. 



(2) 



Then total bending energy of the backbones, Eh = 
J(K/2)[(dt 1 /ds) 2 + (dt 2 /ds) 2 ]ds ||, can be rewritten, 
with the help of Eqs. M and (0), as 
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here L is the total contour length of each backbone. This 
expression proves to be very useful. The second and the 
third terms in Eq.(^|) is deformation energy caused by 
folding of the backbones with respect to the central axis, 
and the first term, n(dt/ds) 2 , is the bending energy of 
DNA central axis contributed by the backbone bending 
rigidity n. So far, basepairs are viewed as thin rods and 
their contribution to the bending rigidity of DNA chain 
is not considered. Because of steric effects caused by 
finite volume and area, basepairs will certainly increase 
the bending rigidity of DNA chain J23| . The simplest way 
to consider such effects is to replace k in the first term 
of Eq. (^) with a phenomenological parameter k* , with 
k* > k. Hereafter this is assumed. 

Besides steric effects, nucleotide basepairs contribute 
also basestacking energy. This energy mainly originates 
from noncovalent van der Waals interactions between ad- 
jacent basepairs |f2l| . Basestacking interaction is short- 
ranged and is characterized by an attraction potential 
proportional to 1/r 6 and a strong repulsion potential 
proportional to 1/r 12 (here r is the axial distance be- 
tween adjacent basepairs). In our continuous model, the 
line density of such Lennare-Jones type potential can be 
written as 



j.[t^£<eo\i2_2(sss<eo.)6] for Up > 0), 

ro LV cos if ' v cos if ' J vr — n 

± [cos 12 ip a ~2 cos 6 (po\ for (cp < 0) , 



(4) 



and the total basestacking energy is Elj = J pds. In 
Eq. (§, r is the backbone arclength between adjacent 
bases; ipo is a parameter related to the equilibrium dis- 
tance between a DNA dimer; e is the basestacking inten- 
sity which is generally base-sequence specific. Here we 
focus on macroscopic properties of DNA and just con- 
sider e in the average sense and take it as a constant, 
with e ~ 14.0fcsT as indicated by quantum chemical cal- 
culations p5]. The asymmetric basestacking potential 
Eq. (|j) ensures a relaxed DNA to take on a right-handed 
double-helix configuration with its folding angle ip ~ ipo. 
However, if adjacent basepairs are pulled apart slightly 
from the equilibrium distance by external forces or ther- 
mal stretching fluctuations, the basestacking interaction 
intensity quickly decreases because of its short-range na- 
ture. In other words, the basestacking potential can en- 
dure only a limited pulling force. We believe this to be 



closely related to the observed DNA highly cooperative 
extensibility. It may also account for the novel elasticity 
of negatively supercoiled dsDNA, since negative super- 
coiling actually leads to an effective pulling force. This 
insight, which is developed in more detail in the follow- 
ing, seems to be confirmed by experiments [ p4| . 

We first discuss the elastic response of the model DNA 
when a pulling force F = /z along direction z is ap- 
plied at its end. The total energy functional is then 
E = Eh+ELj—f^ f cos ipt-zods. And the Green function 
G(t, cp; t', ip'; s) ||,|25|], which determines the probability 
distribution of t and ip. along DNA chain, is governed by 
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where £* = n*/k B T and V(<p) = p(^)/k B T + 
£ p sin ip/R 2 . The spectrum of the above Green equa- 
tion is discrete and hence for long chains, the average 
extension can be obtained either by differentiation of the 
ground-state eigenvalue, g 7 of Eq. (||) with respect to /: 

(Z)/L = (1/L) [ (cosipt ■ z )ds = k B Tdg/df, (6) 
Jo 

or by a direct integration with the normalized ground- 
state eigenfunction, 4>(t,y>), of Eq. (|^): 
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FIG. 1. Force-extension relation of DNA. Experimen- 
tal data is from Fig. 2A of [3]. Theoretical curve is ob- 
tained by the following considerations: (i) £ p — 1.5 nm [4] 



and e = 14.0fc s T [21]; (ii) £* = 53.0/2(cos ip) f=0 



[30], 



r = 0.34/{cos <^}/=o nm and R = (0.34 x 10.5/27r)(tan <p)f=o 
nm [31]; (iii) adjust the value of ipo to fit the data. For each 
ipo, the value of (cosy>)/=o is obtained self-consistently. The 
present curve is drawn with ipo = 62.0° (in close consistence 
with the structural property of DNA [21]), and (cos (p)f=o is 
determined to be 0.573840. DNA extension is scaled with its 
B-form contour length L(cos t^p) /=o ■ 
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Both g and $(t, y>) can be obtained numerically 
through standard diagonalization methods and identical 
results are obtained by Eqs. @ and (Q). The result- 
ing force vs extension relation in the whole relevant force 
range is shown in Fig. [j] and Fig. |^. Our theoretical 
curves are obtained with just one adjustable parameter 
(see caption of Fig. [j]), the agreement with experiments 
is strikingly excellent. According to our theory, the on- 
set of cooperative extension of DNA axial length at forces 
about 70 pN is mainly caused by the yielding of the short- 
range basestacking interaction [^6| . 

Below the onset of cooperative elongation, DNA seems 
to be very stiff and calculations show that at / = 50 
pN the total extension of DNA is only 4.1% longer than 
its B-form contour length, in close accordance with the 
value of 4.6% reported by Smith et al. 0. This is re- 
lated to the fact that the basestacking intensity e is very 
strong p^j . At low forces (/ < 10 pN), because the fluc- 
tuation of the folding angle ip is extremely small, it can 
just be neglected and DNA elasticity is caused by ther- 
mal fluctuations of the axial direction t (entropic elastic- 
ity). It is easy to prove that the now well-known entropic 
elasticity (wormlike chain) model |8| with contour length 
L(cos(p)f—o and persistence length 2£*(cos f) f=o is just 
an excellent approximation of the present theory (here 
(cos(p)/ = o is the average of cosip at zero force). This 
point is demonstrated clearly in Fig. 0. 
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FIG. 2. Low-force elastical behavior of DNA. Here ex- 
perimental data is from Fig. 5B of [1], the dotted curve is 
obtained for a wormlike chain with BPL 53.0 nm and the 
parameters for the slid curve are the same as those in Fig. [l] 

Now we continue to study the supercoiling property of 
the model DNA. Mathematically, a supercoiled dsDNA 
is characterized by its fixed value of linking number Lk. 
It measures the total topological turns one DNA back- 
bone winds around the other or around the central axis, 
and can be expressed as the sum of the twisting num- 
ber, Tiy(ri,r), of backbone ri around the central axis 
r and the writhering number Wr(r) of the central axis, 
i.e., Lk = Tw + Wr According to Eq. (Q), 

Tw(ti,t) = /tx (-b) • d(-b) = J smtpds/R H|7). 
For a linear DNA chain, the writhering number of its 



central axis can be expressed as j28| 

vtr i\ f z x t ■ d(z +t)/ds 

Wr(r) = / — da. 8 

J 1 + z • t 

The elasticity of such a supercoiled DNA chain is de- 
termined by the following energy functional: E — Ef, + 
Elj — f J cos ipt ■ z$ds — TksTLk, where TksT is torque 
associated with the topological constraint. However, the 
writhering number expression given by Eq.(|^) is correct 
only for t • zq ^ — 1, i.e., for chains whose tangential vec- 
tor t never points to — zq [p8| . This condition is satisfied 
actually only for a highly extended chain whose t fluc- 
tuates slightly around Zq. In this case Eq. (|§|) leads to 
Wr(r) ~ (1/2) J(t x dt y /ds — t y dt x /ds)ds, where t x and 
t y are respectively the x and y component of t. This ap- 
proximation is used hereafter. If we are to use Eq.(||) in 
the general case, a cutoff procedure seems necessary to 
avoid divergent results 

The Green equation for this case is written as 

^ - f 92 + 92 + /££ii£t ■ zn - V(w) 
ds — Ue*at 2 + « p cV 2 + k B T 1 z o v (f) 

and the force-extension and torque-supercoiling can then 
be determined through the ground-state eigen rela- 
tionsvalue and eigenfunction of Eq. (||) . Finally, the re- 
lation between extension and linking number is obtained 
by elimination of torque T from these two relations. 

The numerically calculated relations between exten- 
sion and supercoiling degree a at various fixed forces 
are shown in Fig. ^ and compared with the experiment 
of Strick et al. @. Here a is defined by a = ((Lk) — 

(Lk) r=a )/(Lk} r= Q, where (Lfc) r=0 = Jq ds (sin (f) r =a/R 
is the linking number for a torsionally relaxed DNA. The 
parameters for the theoretical curves in Fig. [| are the 
same as those of Figs, [l] and ||, no adjustment has ever 
been made to fit the data. For negatively supercoiled 
DNA, the theory is in quantitative accordance with ex- 
periment (left half of Fig. ||). 
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FIG. 3. Extension vs supercoiling relations at various 
stretching conditions. The parameters for the theoretical 
curves are the same as in Fig. |l| and experimental data is 
from Fig. 3 of [5]. 
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For a < 0, both theory and experiment give three dis- 
tinct regions of DNA elasticity: (i) For forces > 1.3 pN, 
DNA extension does not shrink with the increase of neg- 
ative supercoiling, on the contrary, it may even slightly 
increase as |ct| increases, (ii) For 1.3 > / > 0.3 pN, there 
exists a critical negative supercoiling degree a c . Exten- 
sion of DNA shrinks as a decreases from to a c , then it 
remains approximately constant as a further decreases. 
cr c ~ -0.02 at 0.6 pN. (hi) For / < 0.3 pN, DNA ex- 
tension shrinks constantly with the increase of \a\. In 
this case, no evident difference between the behaviors of 
negatively and positively supercoiled DNAs is observed, 
i.e., DNA can be regarded as achiral 

Thus, the complex elastical property of a negatively su- 
percoiled DNA as well as that of an overstretched DNA 
can be satisfactorily understood by the same framework. 
In this context, although DNA double- helix is quite good 
at enduring external forces it is much weaker at enduring 
torques: while a force ~ 70 pN is needed for a torsionally 
relaxed DNA to trigger cooperative changes of configura- 
tion UK, 0.6 pN is just sufficient for a negatively super- 
coiled DNA with a as small as —2%. This "shortcoming" 
of DNA might have been well noticed and captured by 
various proteins. For example, it seems that RecA pro- 
tein stretches DNA by exerting a torque on the molecule 
& 

However, as shown in the right half of Fig. ^|, for pos- 
itively supercoiled DNA the agreement between theory 
and experiment is poor. It is possible that positive super- 
coiling leads to strong radial as well as axial compressions 
on DNA basepair planes as to make them shrink consid- 
erably or even corrupt. In support of this point, recent 
experiment of Allemand et al. j7| indicates that positively 
supercoiled DNA can take on very surprising configura- 
tions with exposed bases. Therefore, it seems necessary 
for us to take into account the possible deformability of 
DNA basepairs in our theory to understand the elasticity 
of positively supercoiled DNA. We plan to perform such 
an effort. 
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